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Abstract 

The giant halos predicted in neutron-rich Zr isotopes with A = 124 — 138 are investigated 
by using the self-consistent continuum Skyrme Hartree-Fock-Bogoliubov approach, in which the 
asymptotic behavior of continuum quasiparticle states is properly treated by the Green's function 
method. We study in detail the neutron pair correlation involved in the giant halo by analyzing the 
asymptotic exponential tail of the neutron pair condensate (pair density) in addition to that of the 
neutron particle density. The neutron quasiparticle spectra associated with these giant halo nuclei 
are examined. It is found that the asymptotic exponential tail of the neutron pair condensate is 
dominated by non-resonant continuum quasiparticle states corresponding to the scattering states 
with low asymptotic kinetic energy. This is in contrast to the asymptotic tail of the neutron 
density, whose main contributions arise from the resonant quasiparticle states corresponding to 
the weakly-bound single-particle orbits and resonance orbits in the Hartree-Fock potential. 
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I. INTRODUCTION 



The pairing properties in weakly bound nuclei have drawn a lot of attention since the first 
halo phenomenon was discovered in n Li However, the halos observed so far have only 
one or two nucleons in the light exotic nuclei. In order to study the influence of correlations 
and many-body effects, it would be very interesting to investigate the nuclei with a larger 
number of neutrons distributed in the halo. For this purpose, probable halo phenomena 
have been searched for in heavier neutron-rich nuclei Q|. 



The Hartree-Fock-Bogoliubov (HFB) theory [2|-|9j is a powerful tool to describe the heav- 
ier neutron-rich nuclei. It can provide a unified and self-consistent description of both the 
mean field and pairing correlations in terms of the Bogoliubov quasiparticles. Thus, the 
pairing properties in nuclei near the drip-line have been studied extensively within the rel- 
ativistic Hartree-(Fock)-Bog;oliubov scheme (2), I7U18] as well as the non-relativistic HFB 



scheme [4|, [5|, 



13 



19 



-|29|. 



The giant neutron halo with more than two neutrons is predicted first by the relativistic 
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continuum Hartree-Bogoliubov (RCHB) theory IOJ for neutron-rich Zr [2|,[11| and Ca 
isotopes. It is found that the radius of the neutron density distribution shows an abnormal 
increase for A > 60 in Ca and A > 122 in Zr isotopes. It has been shown later that the non- 
relativistic Skyrme HFB model can also describe the giant neutron halo in these elements 



12| and a 



19 



so suppress the asymptotic 



as far as appropriate parameter sets are chosen [13j, [26J. Recently, the giant halo is also 
predicted in Ce isotopes by the relativistic Hartree-Fock-Bogoliubov (RHFB) theory 30]. 

In the preceding works, focuses are often put on the extended tail of the neutron densities 
which leads to the abnormal increase of the root-mean-square (r.m.s.) radius of a halo 
nucleus. It has been discussed that the pair correlation can produce the halo tail of the 
particle density via the continuum coupling 
particle density distribution due to the additional binding 
hand, properties of the pair density or the pairing tensor, which represents the condensate 
of the nucleon pair, especially those in the low-density halo region, have not been studied 
in detail. In this paper, we would like to investigate how the neutron pair density behaves 
in the asymptotic halo region, and also what mechanisms govern its asymptotic behaviors. 

Useful information in analyzing the above properties is the spectrum of the neutron 
quasiparticle states. Because of the shallow Fermi energy, most of the quasiparticle states 



25 



26 



29|. On the other 



2 



are embedded in the continuum energy region by a coupling to the scattering states via the 



pair potential |3j, |4], |6|, |31| . In the standard HFB calculations |2H20 



25 



29], a finite box or a 



harmonic oscillator/Woods-Saxon basis is adopted. In this case the continuum quasiparticle 
states are all discretized (hereafter referred to as the discretized HFB approach), making it 
difficult to describe the asymptotic behavior and to quantify the spectra in terms of, e.g., 



the resonance energy and the width. 



instance, using a very large box 
coupling constant method 
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'here exist techniques to overcome this problem, for 



32 



251 ] . or adopting the analytica 



approach 



33| . or the stabilization method 
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continuation in the 



351 ]. or Gamow HFB 



361 ] . In the present work, we adopt a new formulation of the Skyrme HFB model 



in which the Green's function method 



3jJ is adopted to describe precisely the asymptotic 



behavior of scattering waves for the unbound quasiparticle states in the continuum (hereafter 
referred to as the continuum HFB approach) 27|, |28 ] . 

In this paper we will apply this self-consistent continuum Skyrme HFB approach with 



Green's function method [28J to study the pairing properties in the giant halo Zr nuclei. In 
section II, we will briefly describe the formulation of the continuum Skyrme HFB theory 
with the Green's function method and the numerical details. We will also examine the 
significance of the continuum in the description of the pairing properties by comparing with 
the results obtained by the discretized HFB approach. After presenting the results and 
discussions in Sections III-V we draw conclusions in Section VI. 



II. FORMALISM 

The fundamental building block of the the Hartree-Fock-Bogoliubov (HFB) theory is the 
quasiparticle states. The energy E and the wave function <fri(ra) of a quasiparticle state 
obeys the HFB equation 

( h(ra, r'a') - X5(r - r')5 aa , ~h(ra, r'a') \ , 

dr} A ~ &(rV) = Ei(f)i(ra) 

j \ h*(ra, r'a') -h*{ra, r'a') + X5(r - r')5 atT > I 

(1) 

where A is the Fermi energy. The Hartree-Fock (HF) Hamiltonian h and the pair Hamil- 
tonian h can be obtained by the variation of the total energy functional with respect 
to the particle density matrix p(ra,r'a') = (ip' (r'a')ip(ra)) and pair density matrix 
p(ra,r'a') = {if}(r'a')ip(ra)), respectively. 
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For the spherical system, the quasiparticle wave function can be written as 



(ra 



{r)Yi jm (ra), where 0^(r) 




(2) 



The local particle density p(r) = ^ CT p(ra, ra) and the pair density p(r) = Y2 a p{ r(J , rcr ) 



consist of the products of the quasiparticle wave functions summed up over a 
Using the Green's function technique, the local densities are expressed 



27 



28 



1 the states. 



3 If ] as 



p{r) 
p(r) 



Ij JC E 



(3a) 
(3b) 



where G^fr, r, E) and ^^-(r, r, E) are the radial HFB Green's functions. The Green's 
functions are constructed from the independent solutions of the radial HFB equation ([I]) with 
proper boundary conditions for the wave function (pij(r,E) of the quasiparticle states with 
energy E. Similarly, one can express other local densities needed in the Skyrme functional, 
such as kinetic energy density, spin-orbit densities, etc., in terms of the Green's functions. 
The integrals in Eqs. ([3]) are contour integrals in the complex E plane, and the integration 
path Ce is chosen to be a rectangle with a height 7 = 0.1 MeV and a length E cut = 70 MeV, 
which symmetrically encloses the real negative quasiparticle energy axis as in Ref. 28]. 
The energy step of the contour integration is AE = 0.01 MeV. We choose the box size 
-Rbox = 20 fm, and the mesh size Ar = 0.2 fm for the Runge-Kutta algorithm to obtain 
the independent solution of the radial HFB equation. The quasiparticle wave functions are 



^(out) 



(r,E)/r = (M[ +) (k + r),Bh\ + \k-r] 



connected at r = -Rbox to the asymptotic wave 

where h\ + \z) is the spherical Hankel function 37] and k±(E) = ^2m{\ ± E)/h 2 . Other 
details can be found in Ref. 12811 . 



We employ the Skyrme functional SkI4 



following Ref. 26], with which the giant 

n 

halo phenomenon in the Zr isotopes can be reproduced as predicted in the RCHB theory (2[. 
For the pairing interaction, we adopt the density- dependent delta interaction (DDDI). The 
difference between the zero range and finite range pairing interaction for exotic neutron rich 



nucleus has been discussed in Ref. 

AW = \Vo 



391 ] . The pair field here is taken as 
Pq(r) 



Po 



p(r), q = n or p. 



(4) 
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The parameters in Eg. (U) are adopted as Vq = —458.4 MeV fm 



V 



0.71, a = 0.59, 



and po = 0.08 fm -3 40l443| . which reproduce the experimental neutron pairing gap along 
the Sn isotopic chain. In particular, the parameter Vq is chosen in such a way that the 
DDDI reproduces the 1 S scattering length a = —18.5 fm of the bare nuclear force in the 
low density limit p(r) — > 0, i.e., in the free space outside the nucleus. The quasiparticle 
states are truncated at the maximal angular momentum j max = 15/2 and at the maximal 



quasiparticle energy E cnt = 70 MeV 



26]. 



For the sake of comparison, we also perform a box-discretized HFB calculation. In this 
case the HFB equation ([T]) is solved with the box boundary condition, i.e. by assuming 
that the wave functions of the quasiparticles vanish at the box boundary r = i?box- The 
discretized quasiparticle wave functions thus obtained, <f) n ij(ra), are summed up to construct 
the densities as 



p(r) 
p(r) 



^$>- + l)$>^.(r), 

Ij n 



(5a) 
(5b) 



III. PAIR CORRELATION IN GIANT NEUTRON HALO NUCLEUS 

The filled circles in Fig. []7a) are the two- neutron separation energy S2 n {N,Z) = 
E(N,Z) — E(N — 2,Z), where E(N,Z) is the total binding energy of the isotope with 
N neutrons and Z protons obtained in the continuum HFB calculation. It suddenly drops 
from S% n = 4.81 MeV at 122 Zr to 5*2 n = 0.64 MeV at 124 Zr as the neutron number exceeds the 
magic number N = 82. Then it gradually decreases to an extremely small value 0.04 MeV 
at 138 Zr, and finally becomes negative in 140 Zr. Thus 138 Zr is the neutron drip-line nucleus 
in the present model, which is consistent with the previous investigation 26j. The small 
two-neutron separation energy is one of the conditions for the emergence of the halo struc- 
ture. Figure [TJ (b) shows the neutron Fermi energy A as well as the Hartree-Fock (HF) 
single-particle energies e, which are the eigen energies of HF Hamiltonian h (obtained after 
the final convergence of the continuum HFB calculation). One can see that, as the neutron 
number increases, the Fermi energy is raised up, to a position quite close to the continuum 
threshold, while all the HF single-particle orbits fall down. Explicitly, the 3pi/2, 3p 3 / 2 and 



2/7/2 states evolve from unbound resonances (A < 122) to weakly bound orbits (A > 126), 
and the 2/5/2 orbit remains as a resonance in the continuum. The values of these HF single- 
particle energies in 124 Zr ~ 138 Zr can be found in Table [H The HF Hamiltonian has another 
resonance lh 9 / 2 which lies slightly higher than the plotted energy region in Fig. [T] (b). 

Figure H] (a) shows the neutron density p(r) calculated for the Zr isotopes with A = 122 — 
138 using the continuum HFB calculation (solid lines). It is seen that the neutron density 
for A = 124 — 138 exhibits a long tail extending far outside the nuclear surface compared 

1 /2 

with A = 122. Figure [2] (b) shows the r.m.s. radius r r m . s . = [J 47rr 4 p(r)dr / J 47rr 2 p(r)c?r] 
calculated for the corresponding neutron density (filled circles). Compared with the isotopic 
trend in A < 122, which gives an extrapolation as r r m Si ps 0.87N 1 / 3 fm, the neutron r.m.s. 
radius in 124 Zr and the heavier isotopes displays a steep increase with N. These results 
are consistent with the previous investigations in Refs. [2| and 26|, where the giant halo 
structure, i.e. the long tail which can accommodate more than two neutrons, is predicted 
in the isotopes with A > 124. 

Figure [2](c) shows the neutron pair density p(r). We clearly see that the pair density has 
an even more significant tail, which is characterized by i) p(r) > p[r) for r ^ 10 fm, and ii) 
that p(r) has a more gentle slope than p(r). These features are reflected in the r.m.s. radius 

1 /2 

weighted with the neutron pair density, f r m . s . = [J Anr 4 p{r)dr / J 4iir 2 p(r)dr^ , which is 

plotted in Fig. [2] (d). We shall call r r m s the pair r.m.s. radius hereafter. It is significantly 

larger than r r . m . s . of the neutron density. Also, the neutron pair r.m.s. radius f r . m . s . exhibits 

a sudden jump at 124 Zr, and remains almost the same until the drip-line. Here, f r . m . s . of 

122 Zr is omitted in the plot due to the absence of pairing. 

It is argued in Refs. and 26] that the neutron orbits 3pi/ 2 , 3p 3 / 2 , 2/ 5 / 2 and 2/ 7 / 2 

around the Fermi energy play a decisive role in forming the giant halo structure in the 

neutron density. Similarly to Refs. Q and 26], we decompose the neutron density with 

respect to different partial waves Ij, and calculate separately the r.m.s. radius 

f Anr 4 pij{r)dr\ 1 ^ 2 
J Anr 2 pij{r)dr J 

weighted by the corresponding /j-decomposed neutron density 

(2j 



r 



r.m.s. ,lj 



(6) 



Pij{ r ) 



which sums the contributions from the quasiparticles with the quantum number Ij within 
E = ~ E cut . The r.m.s. radii for Ij = s 1/2 , P1/2, P3/2, 4/2, d 5/2 , / 5 / 2 and / 7/2 partial 
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waves are plotted (filled symbols) in Fig. |3j We can see clearly that the increase of the total 
r.m.s. radius is mainly contributed by the p and / states, in agreement with Refs. 0, [26^ 

The pair density displays an apparently different behavior from that of the neutron density 
as seen in Fig. HJ in which we plot the neutron pair r.m.s. radius 

J 4-irr 4 pij(r)dr\ 



r.m.s.,Zj 



J 4:7rr 2 pij(r)dr 
weighted by the /j-decomposed neutron pair density 

(2j 



(8) 



Pij{ r ) 



2 i^^i dES ^r,r, E) . (9) 

A notable feature is that the jump at N — 124 is seen not only for the partial waves 
Pi/2,P3/2, fb/2 and /V/2, which have HF single-particle orbits or resonances near the Fermi 
energy, but also for the partial waves £1/2,^3/2 and This indicates that the s- and 

<i-waves also play important roles in the large tail of the pair density. 

To clarify the difference, we present in Fig. [5] the composition of the neutron density p(r) 
and pair density p(r) with respect to the partial waves I J, i.e., py(r)/p(r) and Pij(r)/p(r) for 
126 Zr and 138 Zr. As shown in the previous works 2[ and 26], the most dominant compositions 
of the neutron density p(r) in the halo region ri^lO fm are the p-waves, and the next 
are the /-waves, while the s- and d-waves as well as others partial waves have very little 
contributions, and almost vanish at large distances r ^ 15 fm. 

For the pair density p(r), however, we can see from Fig. 5 (c) and (d) that the partial 
wave composition is very different from that of the neutron density p(r). The dominance 
of the waves P1/2, P3/2, /s/2 an d /V/2 is commonly seen also for the pair density, but to a 
less extent. The waves £1/2,^3/2 and 0^5/2, on the other hand, have small but non-negligible 
contributions in the halo region r ^ 10 fm. This feature becomes even stronger in 138 Zr than 
in 126 Zr. Since there is no bound HF single-particle orbits nor resonance orbits in the s- and 
(i-waves near the Fermi energy, it is non-resonant continuum orbits that contribute to the 
pair density. 

Before closing this section, we would like to compare the present results obtained in the 
continuum HFB calculation with those obtained in the box-discretized HFB calculation. In 
the latter case, the quasiparticle states are discretized even in the continuum energy region 
E > |A|, and the densities are obtained as Eq. flH]). 

The two-neutron separation energy S^™ in this approximation is shown in Fig. [1] (a) with 
the open circles. The total neutron density and pair densities of the box-discretized results 



are shown with the dotted lines in Fig. [2] (a) and (c). The corresponding neutron r.m.s. 
radii r r . m . s . and pair r.m.s. radii f r . m . s . are shown with open circles in panels (b) and (d). 
The r.m.s. radii of different partial waves are also shown in Figs. Eland H] with open symbols 
for the neutron density and the pair density, respectively. 

From the comparison between the continuum and the box-discretized calculations, we 
observe the following. The influence of the box boundary condition is seen both in the 
neutron density p(r) and in the pair density p(r) at r « 15 — 20 fm close to the box boundary 
(Fig. [2] (a) and (c)). However this influence causes only negligibly small difference in the bulk 
properties such as the total energy and its derivative, the two-neutron separation energy S^n. 
This is because the density and the pair density of neutrons at r ^ 15 fm are very small, and 
its contribution to the total energy is accordingly small. The r.m.s. radius r r . m . s . (Fig. [21(b)) 
of the total neutron density is also affected very little by the box-discretization. If we look 
into the r.m.s. radii r r m s ^ for each partial waves (Fig. [3]), we can see some difference in the 
contributions from the p-wave while the other partial waves are not affected. This indicates 
that the neutron density contributed from the p-wave is the most extended while other 
partial waves extend less (see also Figj5]). More importantly, we can see more significant 
difference in the neutron pair r.m.s. radii r r m s and r T ms ^j (cf. Fig. [2] (d) and Fig. [4]). 
The difference is seen not only in the p- and /-waves, but also in the s- and <i-waves. This 
difference is larger for lower angular momentum states and in more neutron-rich nuclei, 
where we expect more continuum coupling due to the lower centrifugal barrier. It can be 
concluded from these observations that the proper treatment of the continuum quasiparticle 
states is important for the description of the tail part of the pair density, i.e. the pair 
correlation in the giant halo. 



IV. CONTINUUM QUASIPARTICLE SPECTRA 



A. Low-lying quasiparticle states 



A quasiparticle state with excitation energy E in the partial wave Ij has contributions 
to the density p(r) and the pair density p(r), which are given by 27 . 



Airr 2 



r, r, 



7T 



(2j + 1) 1 
4irr 2 ir 



Im£ ( 



o,ij 



-E 
-E 



(10) 



S 



with the use of the HFB Green's function and an infinitesimal constant e. They satisfy 
p(r) = I dEpij(r, E) and p(r) = I dEpij(r, E). We can then define the occupation 
number density as a function of E by 



which represents the neutron number associated with the quasiparticle state with energy E 
and the quantum numbers Ij. We can also define the pair number density 



using the pair density contribution (r, E) . 

Figure E] (a)-(d) show the occupation number density riij(E) of neutron quasiparticle 
states in a quasiparticle energy interval E = ~ 2 MeV for 126 Zr and 138 Zr, while Fig. O 
(e)-(h) are the neutron pair number density nij{E) for the same isotopes. We choose e = 
10~ 10 MeV in Eq. (fTUj) . Note that the quasiparticle states with energy E larger than the 
threshold energy |A| form a continuum spectrum. 

We can see from this figure that the quasiparticle spectra of 126 Zr and those of 138 Zr 
have common features. The partial waves Pi/2,P3/2, h/2 and fan in panels (c), (d), (g) 
and (h) have peak structures with finite width, i.e., quasiparticle resonances, in the spectra 
of nij(E) and hij(E). These quasiparticle resonances correspond to the HF single-particle 
orbits which are shown in FigJT] (b). It should be noticed here that in 126 Zr the HF single- 
particle orbits 3^3/2 and 2/7/2 are discrete bound states with very small binding energy 
e ~ —0.1 ~ —0.2 MeV. It is the pair correlation that transforms these bound HF single- 
particle orbits to unbound quasiparticle resonances. Similarly, the weakly bound HF single- 
particle orbits 3p3/2,3pi/2 and 2/7/2 in 138 Zr become the quasiparticle resonances when the 
pair correlation is taken into account. Table 1 lists the resonance energies E res of these 
quasiparticle resonances, which we evaluate as the peak energy of nij(E). 

Another noticeable feature is seen in the quasiparticle spectra of the partial waves with 
positive parity, Si/ 2 , d 3 / 2 , e/5/2, etc, plotted in Fig. 0(a), (b), (e) and (f). We see here smooth 
profiles of nij(E) and nij(E), indicating non- resonant continuum quasiparticle states. Con- 
cerning the occupation number density riij(E), the positive parity non-resonant quasiparticle 
states have contributions which are smaller, by several orders of magnitude, than the con- 
tributions of the resonant quasiparticle states with the negative parity shown in panels (c) 




(11) 




(12) 
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and (d). On the other hand, the pair number densities nij(E) of the positive parity partial 
waves ((e) and (f)) are more than ten times larger than nij(E) of the same partial waves 
((a) and (b)). It is also noticed that nij(E) of the positive parity partial waves ((e) and 
(f)) are even comparable to those of the negative parity partial waves ((g) and (h)), except 
around the resonant peaks. These sizable contributions from the non- resonant S1/2, d%/2 and 
^5/2 quasiparticles to the pair number density can be related to the non- negligible fraction 
pij{r)/p{r) of the £1/2,^3/2 and d 5 / 2 quasiparticles in the halo tail of the pair density as 
shown in Fig. [5] (c) and (d). 



B. Pairing gaps and missing bound quasiparticle states 

In the isotopes 126_138 Zr, all the quasiparticle states are embedded in the continuum 
energy region E > |A| (Fig. [6] is an example), namely there exist no bound quasiparticle 
states with E < |A| except in 124 Zr. [In 124 Zr there is only one bound quasiparticle state 
3p 3 /2 at E = 0.436 MeV, which is very close to |A| = 0.446 MeV.] It is a remarkable 
effect of the pair correlation since the HF single-particle states 3pi/2,3p3/2 and 2/7/2 are all 
bound orbits if the pair correlation were neglected (cf. Fig. QJb)). To understand this, we 
first evaluate the average pairing gaps defined as usual by A uv = J drA(r)p(r) / J drp(r) 
and A vv = J drA(r)p(r)/ J drp(r) with two different weight factors given by the neutron 
density or the neutron pair density. They are listed in TableHl Their values A uv m 0.59 — 0.68 
MeV and A vv « 0.53 - 0.61 MeV are fairly constant (A « 0.5 - 0.7 MeV) in the interval 
A = 126-138, except for A « 0.42-0.47 MeV at A = 124, where basically only two valence 
neutrons participate in the pair correlation. We observe here that the average pairing gap 
A ?a 0.5 — 0.7 MeV is larger than the absolute value of the Fermi energy |AJ_in the isotopes 
126 ~ 138 Zr, where |A| < 0.42 MeV. This relation A > |A| is proposed in Ref. 4j as a criterion 
for a non-perturbative role of the pairing on the properties of the weakly bound nuclei. We 
can also examine this criterion with respect to the individual quasiparticle resonances. We 
here note that the energy of the quasiparticle resonances 3pi/2 and 3^3/2 satisfy E TCS > 0.4 
MeV and E rDS > |A| in the isotopes with A = 126 — 138 (and similarly, E rcs > 0.7 MeV and 
Ercs > |A| for 2/7/2). As a result, if we use the BCS expression E ms w — A) 2 + A 2 ff 
to relate the quasiparticle energy E rcs and the effective pairing gap A e g relevant to the 
specific quasiparticle state, we can estimate the effective pairing gap of the 3pi/ 2 and 3p 3 / 2 



10 



resonances as A c ff <; 0.4 MeV (and A c g <: 0.7 MeV for the 2/7/2 resonance) and A e g- ^ |A|. [ 
Here we used also the fact that the HF single-particle energies of 3pi/ 2 ,3p 3 / 2 and 2/7/2 are 
almost degenerate with the Fermi energy, e ~ A (cf. Fig.[Tj).] Again we can infer the relation 
A e ff > |A| also for the effective pairing gap of the weakly-bound orbits under consideration. 



C. Width of quasiparticle resonances 

We evaluate the width of a quasiparticle resonance by reading the FWHM of the peak 
structure in the occupation number density riij(E). The resonance energies and widths of 
the low-lying quasiparticle resonances 3pi/2, 3p3/2, 2/5/2 and 2/7/2 are displayed in Fig. [7J 
The length of the bar in this plot represents the width T of the resonant states multiplied by 
a factor of 5. For the sake of reference, it also plots the threshold energy |A| and the position 
Vb of the barrier top of the HF potential including the centrifugal potential measured from 
the Fermi energy. We list in Table 1 the barrier height of the potential V max = Vb — |A| 
and the resonance energy e res = E res — |A| = E TCS + A measured from zero HF single-particle 
energy. 

There are at least two mechanisms that govern the width of the quasiparticle resonances. 
One is the barrier penetration of a single-particle motion, which is present in any potential 
models, i.e. the HF potential in the present case. The other is the one caused by the pair 
correlation, through which even bound HF single-particle orbits can couple to the scattering 
states. An ex amp le of the latter is known as quasiparticle resonances originating from deep 
hole states j^, 31]. 

These two mechanisms are interwoven in a rather complex way to produce the widths 
of the low-lying quasiparticle resonances under discussion. A typical example is seen in the 
isotopic dependence of the width of the 3pi/2 resonance (Fig UJ (a)). As the mass number 
increases from A = 124 to A = 138, the width of the 3pi/2 resonance first decreases, then 
turns to increase at larger neutron number A ft 132. The decrease of the width in the interval 
124 < A^ 132 is hardly explained by the barrier penetration mechanism alone: the barrier 
height Vmax decreases (see Table [I]) while the resonance energy e rcs are fairly constant in this 
interval, thus the barrier penetration mechanism would lead to an "increase" of the width. 
The decrease of the width may be explained by combining the potential barrier penetration 
and the continuum coupling caused by the pairing in the following way. We here note that 
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because of the change in the relative ordering of the HF single-particle energy e and the 
Fermi energy A (cf. Fig.2), the particle-character of the 3pi/2 quasiparticle state, dominant at 
A = 124, is weakened gradually with increasing the mass number from A = 124 to A 132, 
and a hole-character grows with further increase of A. If we can assume that the coupling 
of the hole-component y?2,zj( r ) to the scattering wave via the pair potential is weaker than 
that of the particle-component <pi,ij(r) to the scattering wave via the barrier penetration, 
we then expect the decrease of the width. At the mass numbers A « 132 — 134, the barrier 
penetration becomes progressively effective since the resonance energy e res increases, and it 
causes the increase of the width. However, the width remains finite even at A = 136, 138 
where the resonance energy lies above the barrier height. In these isotopes, the width is 
probably controlled dominantly by the coupling to the continuum via the pair potential. 
These speculations remain at a qualitative level. The analytical evaluation of the width of 
the quasiparticle resonance is known only for those originating from deep hole states j^, 31], 
but our situation is more complex. Quantitative understanding of the width of the low-lying 
quasiparticle resonances associated with the weakly-bound single-particle orbits remains as 
a future subject. 



D. Continuum coupling and comparison with the box discretization 

It is interesting to compare the box-discretized HFB calculation and the continuum HFB 
calculation for the description of the resonant quasiparticle states. 

The quasiparticle states are all discretized when the box boundary condition is adopted, 
and the discretized states corresponding to the quasiparticle resonances we are discussing are 
the lowest energy states in each of the partial waves Pi/2,P3/2> h/2 and iV/2- The quasiparticle 
energies of these states are shown in Figf7]with open circles. We see large deviations (the 
deviation is even larger than the value of the resonance width) from the resonance energies 
E rcs obtained in the continuum HFB calculation for the 3pi/2 and 3^3/2 states at A ^134. 
Contrastingly, at A = 136 and 138, the quasiparticle energies obtained in the box-discretized 
calculation agree rather well with the resonance energies (the deviation is comparable with 
the value of the resonance width). A good agreement is also seen in the quasiparticle 
resonances of 2/ 7 / 2 and 2/ 5 / 2 . 

A hint of the large deviation is suggested in the ratio between the width T and the 
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resonance energy e rcs . It is seen from Table [J that the quasiparticle resonance 3pi/2 at 
A = 124, for instance, has a width T = 61 keV which is comparable to the resonance energy 
e re s — 79 keV, i.e., r/e ros ~ 1. This ratio decreases with A. We find that the large deviation 
is seen in the case of the large width-energy ratio r/e res <^ 1/3 (at 124 < A & 134). The same 
is seen also for the 3^3/2 resonance. On the other hand, the small deviation can be linked 
to the small width-energy ratio r/e rcs ^l/5 as seen in the 3p±/2 and 3p^/2 resonances at 
A = 136 and 138, and also in the 2/5/2 and 2/7/2 resonances in all the isotopes. Clearly the 
large width-energy ratio r/e res ^ 1/3 indicates that the coupling to the continuum scattering 
states is strong even though the resonance energy is located below the potential barrier 
(e res < Vmax or E < Vb)- Consequently the bound state approximation does not work. For 
the same reason it is hard to describe these resonances in a box-discretized calculation unless 
the box size is taken sufficiently large. In order to describe the resonance at such a small 
energy e res ~ 50 keV, one has to use a large box satisfying h 2 /2m x (vr/i? box ) 2 ^ e res , i.e., 
-Rbox ~ 60 fm. A further larger box is needed in order to describe the distribution around 
the peak [jj]] to evaluate the width, or one needs to combine with other methods, such as 



stabilization method 



35j, etc. 



From the above analysis, we can conclude that the correct treatment of the continuum 
is important to describe the quasiparticle resonances originating from weakly bound orbits 
with low angular momenta. 



V. EXPONENTIAL TAILS OF THE PAIR CORRELATED HALO 

In this subsection, we will discuss the asymptotic behavior of the density and the pair 
density associated with the giant neutron halo. We have already discussed in section III 
that the density and the pair density exhibit extended tails which can be characterized by 
an exponentially decreasing behavior as a function of r. In order to understand the origins 
of the tail behaviors, we shall investigate the neutron density pij(r) and the neutron pair 
density pij(r) which are decomposed with respect to the partial waves Ij. 

Figure 8 (a) and (b) shows the Zj'-decomposed neutron density 47rr 2 p^(r) weighted with 
the volume element 4irr 2 for 126 Zr. We already saw the dominance of the p-wave component 
in the halo tail region (r ^ 10 fm) in connection with Fig. in section III. Noticeably the 
exponential slopes in the tail region are different for different partial waves. The pi/2 and 
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£>3/2 components have the most gentle slopes (hence they dominates in the tail) while the 
slopes of fj/2 and f§/2 are steeper than the p- waves. The densities arising from the s- and 
<i-waves are much steeper, and thus their contributions to the halo tail is negligible. 

Looking at the /j-decomposed neutron pair densities 47ir 2 pij(r), shown in panels (c) 
and (d), we immediately see that the exponential slopes of the neutron pair densities are 
apparently different from the slopes associated with the neutron density pij(r) (panels (a) 
and (b)), and it is also obvious that the exponential slopes exhibit much weaker dependence 
on Ij. The slopes are rather similar among the partial waves S1/2, (£3/2 and as well as 
Pi/2,P3/2,/5/2 and f 7/2 . 



A. neutron density 

In order to clarify the behaviors of the exponential slopes in the halo tail region, we shall 
quantify the exponential slope of the neutron density pij(r). For this purpose, we consider 
a simple fitting function 

r 2 pRt,ij(r) oc (rh\ + \inij r)\ (13) 

with h\ (z) being the spherical Hankel function 37]. A concrete form for the orbital an- 

r 1 2 

gular momentum / = 1 is r 2 p& t) ij{r) = Ce~ 2Rl]r 1 H — — . This is an asymptotic form 
expected to be obtained if the tail density is contributed by a single quasiparticle state with 
a fixed energy. (Note that the asymptotic form of the second-component wave function is 
(p2,ij(r)/r ~ h\ ; (i«yr) with «y = ^2m(E + \\\)/h 2 .) We fit this function r 2 p^ t ^{r) to the 
numerically obtained density r 2 pij(r) in the tail region, and we extract the parameter 2nij 
that represents the exponential slope of r 2 pij(r). A simpler choice of the fitting function 
would be r 2 pfi t jj(r) = C exp(— 2nijr) valid for n^r ^> 1. But we do not adopt this because 
the effect of the centrifugal barrier is not completely negligible (i.e. «yr 3> 1 is not realized) 
in the region under discussion, r ~ 20 fm. We use the interval r = 15 — 20 fm for the fit, 
and we denote the extracted value as 2k^, which we shall call as asymptotic exponential 
constant. 

The extracted asymptotic exponential constant 2kij is shown in Fig. |9] (a) and (b). It is 
noticed that the extracted values 2kij differ for different Ij, and their isotopic dependencies 
are also different. For instance, 2/cy for / 5 / 2 is obviously different from that for / 7 / 2 although 
both equally have sizable contributions to the halo density. 
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We shall now show that the asymptotic exponential constants of these partial waves 
are governed by the low-lying quasiparticle resonances 3pi/ 2 , 3p 3 / 2 , 2/ 5 / 2 and 2/ 7 / 2 . To show 
this, we evaluate a part of the neutron density pij{r), denoted by p'^ir) below, in which we 
consider only contributions from the low-lying quasiparticles with E = — 3 MeV including 
the resonances mentioned above. The truncated neutron density p[j(r) can be calculated by 
using Eq. (J7J) in which the contour of the integral encloses only the quasiparticle states with 
E = — 3 MeV. In Fig. [10] (a) we compare p^-(r) with pij(r) for p 3 / 2 contributions in 138 Zr. 
It is seen that p' P3/2 { r ) reproduces p p . i/2 (r) very nicely as far as the tail region r > 10 fm is 
concerned. We extract the asymptotic exponential constant 2k^ for p' Pa/2 {f) by using the 
same fitting function of Eq. (113j) . and we confirm that 2Ky and 2k^ are almost the same. 
We also compare the tail of the neutron density pij(r) with that of the density pij(r,E res ) 
which represents a contribution from the quasiparticle state at the resonance energy E res 
(cf. Tabled]). It is seen in Fig. [10] that the exponential slope of pij(r,E Tes ) agrees quite well 
with that of pij{r). 

To be more quantitative, we compare in Fig. [9] (b) the extracted asymptotic exponential 
constant 2k^ of pij(r) with the exponential constant 

/ 2m(ff rcs +|A|y 
2n lj:Ves = 2y — (14) 

of pij(r, -Eres), which can be calculated using the resonance energy E TCS , for the partial waves 
Pi/2)P3/2> fb/2 and f 7 /2 in all the isotopes. The agreement is very good and it is even hard 
to see the difference between the two in some cases. We thus confirm that the low-lying 
quasiparticle resonances 3pi/ 2 , 3p 3 / 2 , 2/ 5 / 2 and 2/ 7 / 2 corresponding to the weakly-bound HF 
single-particle orbits dominate the tail and govern its asymptotic exponential constant. 

The asymptotic exponential constants 2nij of the si/ 2 , cZ 3 / 2 and d§/2 partial wave neutron 
densities behave in a different way as seen in Fig. [9] (a). Firstly, we observe the large values 
of 2nij, which can be attributed to the fact that the tail of pij(r) in these partial waves is 
dominated by deeply bound orbits. Secondly we find that the values of 2k^ decrease steeply 
with increasing A. This feature can be explained by the fact that the non-resonant continuum 
states in the low-lying region of these partial waves contribute more as the neutron drip-line 
is approached. In Fig.lTOl(b) we plot the truncated neutron density p'lAr) of syz partial wave 
in 138 Zr, which represents the contributions from low-lying non- resonant continuum states. 
Although p'lAr) is smaller than pij(r) by several orders, it dominates the far external part 
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(r > 15 fm) of the total density pij(r). Note, however, that the si/ 2 , d 3 / 2 and d 5 / 2 densities 
do not influence the tail of the total neutron density as their relative contributions are 
negligible in comparison with those contributed from the low-lying p, /-wave quasiparticle 
resonances. Therefore, the complex behavior of 2«;y in the s, ci-wave does not affect the 
conclusion on the dominance of the low-lying quasiparticle resonances in the exponential 
tail of the neutron densities. 



B. neutron pair density 

Let us now examine the asymptotic exponential constants of the neutron pair densities 
pij{r). Similarly to the neutron density Py(r), we fit an exponentially decreasing function to 
the numerically obtained neutron pair density py(r), but we need a proper fitting function 
other than pfit,y(r) (Eq. (ll3j) ) which is appropriate only for the normal density. We adopt 

r 2 p&t,ij{r) oc rh\ + \iRijr) (15) 



1 + ^- 



as a fitting function. Specifically, for / = 1, we have r 2 p&t,ij(r) = Ce~ Klj 
This functional form is based on the following considerations. We first point out that 
there is no simple estimate for the asymptotic form of the pair density. Even if we as- 
sume that a quasiparticle state with energy E contributes to the pair density, the con- 
tribution r 2 pij(r,E) oc (pi,ij{r)(f2,ij(r) does not have an exponential form since the first- 
component wave function <£>i,y(r) is an oscillating sinusoidal function in the asymptotic 
region as the quasiparticle states under consideration are all embedded above the thresh- 
old energy. We therefore exploit only the asymptotic form of <P2,ij(r)i which behaves as 
^2,ij( r ) ~ rh[ + \iKijr), to prepare the fitting function. We use the same interval r = 15 — 20 
fm for the fitting, and denote the extracted asymptotic exponential constant as «iy for the 
pair density. 

The extracted value of the asymptotic exponential constant Ky is plotted in Fig. [9] (c) 
and (d). It is seen here that the behaviors of «y are different from those of the asymp- 
totic exponential constants 2«y of the neutron densities Py(r). This immediately demon- 
strates that the quasiparticle resonances are not the major origin of the exponential tail 
of the pair densities Py(r). It is also seen that the values of Ky for all the partial waves 
s i/2,Pi/2,P3/2, d 3 /2, d 5 / 2 , fs/2 and fa/2 are rather similar, and they all exhibit similar isotopic 



1(3 



dependencies. 



It is useful here to refer to the argument in Ref. [4j which asserts that the asymptotic 
exponential constant of the pair density may be given by 



K^=yj-±± ■ (16) 

This is based on the consideration that the asymptotic behavior of the pair density is gov- 
erned by the contribution of the lowest energy quasiparticle state, and in the case of nuclei 
near the drip-line, the lowest energy quasiparticle state is the continuum state with E — |A| 
(cf. Eq. (I14p ). In fact we see in Fig. [9] (c) and (d) that the isotopic dependencies of Rij has 
some resemblance to that of k m \ n although a difference by a factor of about 2 is seen here. 
We thus deduce that non-resonant continuum quasiparticle states dominates the asymptotic 
tail of the pair density,although the difference by a factor of ~ 2 suggests something beyond 
the argument in Ref. |4|. 

Figure ITOl(c) and (d) are made to check the above deduction. In these panels we compare 
the neutron pair density pij(r) with the truncated pair density p^(r) where only the low- lying 
quasiparticle states with E < 3 MeV are included in the sum of Eq. (JTj). It is also compared 
in panel (c) with the pair density pij(r,E res ) that originates from the single quasiparticle 
state at the resonance energy E res for p 3 / 2 partial wave. It is seen that pij(r, E res ) apparently 
fails to describe the exponential tail since <fijj(r, E rcs ) and hence pij(r, E res ) oscillate with r. 
On the other hand the exponential tail of the pair density pij (r) is nicely reproduced by the 
pair density pj -(r) truncated with E < 3 MeV. This observation applies to the pair densities 
of all the p, f and s, d waves (see Fig. [TO] (d) for S1/2). Recall that the s, d- waves have only 
non-resonant quasiparticle states in the low-lying spectrum E < 3 MeV. We thus conclude 
that the non-resonant continuum quasiparticle states dominate the asymptotic exponential 
tail of the neutron pair density. It is in contrast to that of the neutron density in which the 
low-lying quasiparticle resonances dominate. It is also concluded that the exponential tail 
of the pair density is governed by coherent superpositions of the contributions from the non- 
resonant continuum quasiparticle states near the threshold energy. Under this situation, we 
can evaluate approximately the exponential tail arising from the non-resonant continuum 
quasiparticle states as 

r 2 p lj (r)= [ dE<p hlj (r,E)<p 2 , l:i (r } E)~ I desm(kr)e- Kr (17) 
J\x\ Jo 
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with k = ^2me/h 2 1 k = y^2m(2\X\ + e)/H 2 and e = E — |A|. Since the weight factor e~ Kr 
is dominant in an interval < e^|A|, the superposition of non-resonant continuum wave 
sin(kr) with such a weight factor gives ~ e _Knrr with k QT ~ 0(k min ). Consequently we can 
expect pij{r) ~ e~ KljT with Rij = K miQ + k nr ~ 2/€ min . This estimation is in qualitative 
agreement with the observation kij rj 2K m ; n seen in the extracted asymptotic exponential 
constant «y (cf. Fig. [9] (c) and (d)). 

VI. CONCLUSIONS 

We have investigated the pair correlation in 124 - 138 Zr by using the self-consistent contin- 
uum Skyrme Hartree-Fock-Bogoliubov (HFB) approach, in which the asymptotic behaviors 
of the wave functions of the continuum quasiparticle states are properly treated with the 
Green's function method. The giant neutron halos in these very neutron-rich weakly-bound 
nuclei are analyzed in detail. Focuses are put on properties of the halo part of the neutron 
pair condensate, i.e. the exponential tail of the neutron pair density p(r), in particular, 
the r.m.s. radius, its single-particle composition, and the exponential slope parameter (the 
asymptotic exponential constant). We found that these are apparently different from the 
corresponding quantities associated with the exponential tail of the neutron density p(r). 

In order to clarify the origin of the differences, we have looked into the spectrum of 
the neutron quasiparticle states. Because of the small Fermi energy A and of the relative 
largeness of the pairing gap A^|A|, all the quasiparticle states in 126 ~ 138 Zr are located 
above the threshold energy. Typical examples are the quasiparticle resonances 3pi/ 2 , 3p 3 / 2 
and 2/ 7 /2, which play central roles in forming the giant neutron halo. They originate from 
the weakly bound orbits in the Hartree-Fock (HF) potential, but they appear as resonances 
with finite width once the pair correlation is taken into account. There exist non-resonant 
continuum quasiparticle states, and they contribute also to the giant halo. Such quasiparticle 
states are seen in all the partial waves, and typically in S1/2, g?3/2 and d 5 / 2 - 

The central finding in the present analysis is that the exponential tail of the neutron pair 
density pij{r) is contributed mainly from the low- lying non- resonant continuum quasipar- 
ticle states. This is contrasting to the microscopic structure of the exponential tail of the 
neuron density pij{r), the main contributions of which are from the low- lying quasiparticle 
resonances 3pi/2,3p 3 / 2 and 2/ 7 / 2 as well as 2/5/2, all corresponding to the HF orbits and 
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resonances near the Fermi energy. 

The different microscopic origins of the neutron density pij(r) and of the neutron pair 
density pij{r) reflect to the asymptotic exponential constants 2«y and R\j parameterizing 
the exponential slope of the tail of the neutron density pij(r) and the neutron pair density 
Pij{r), respectively. Concerning the neutron density pij(r) of the dominant partial waves px/2 
and P3/2, subdominant /s/2 and ft/2, we found that the constant 2nij is governed essentially 
by the peak energy of the lowest-lying quasiparticle resonance that arises from the bound or 
resonant HF single-particle orbit closest to the Fermi energy. The asymptotic exponential 
constant of the neutron pair density pij{r) has a close relation to the threshold energy |A| 
of the non-resonant continuum states, rather than the energies of the low-lying quasiparti- 
cle resonances, since the superpositions of the non-resonant continuum quasiparticle states 
dominate the tail of pij{r). 

In the present analysis, we have also put a focus on the comparison between the continuum 
HFB calculation with the box-discretized HFB calculation. This comparison has pointed 
to a few specific problems of the box calculation, which arise if the size of the box is not 
sufficiently large. One is that the finite box causes relatively large error in describing the 
exponential halo tail of the neutron pair density pij(r), in comparison with the tail of the 
neutron density pij{r). This is because the halo tail of the neutron pair density pij(r), which 
is dominated by the contributions of the non-resonant continuum quasiparticle states with 
small asymptotic kinetic energies e = E — |A|, is more extended than the neutron density 
pij(r). Another problem is found in describing the quasiparticle resonances corresponding 
to the weakly-bound 3pi/2 and 3£>3/2 orbits, in particular in the isotopes with 132. We 
here have only small centrifugal barriers in the p-waves, which makes the coupling to the 
scattering states strong. Consequently, the quasiparticle resonances 3pi/2 and 3p^/2 are not 
represented well by any eigen states in the discretized calculation with a small box. 
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FIG. 1: (a) Two neutron separation energy S^n obtained in the continuum (filled circle) and box- 
discretized (open circle) Skyrme Hartree-Fock-Bogoliubov (HFB) calculations for neutron-rich Zr 
isotopes with A = 118 — 138. (b) Neutron Hartree-Fock (HF) single-particle energy e of Zr isotopes 
around the Fermi energy, which is given by the eigen solution of HF Hamiltonian h after the final 
convergence of the continuum HFB calculation. The dashed line denotes the Fermi energy A. 



22 




5 10 15 20 116 120 124 128 132 136 140 



r [fm] A 

FIG. 2: (a) Neutron density p(r), (b) neutron root-mean-square (r.m.s.) radius r r . m . s ., (c) neutron 
pair density p(r) and (d) neutron pair r.m.s. radius f r . m . s . for Zr isotopes. The solid lines and the 
filled circles are the results of the continuum HFB calculation while the dotted lines and the open 
circles are those obtained in the box-discretized HFB calculation. 
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FIG. 3: Neutron root-mean-square radius r r . m . s y of the s, p, d and / partial waves of Zr isotopes 
calculated for the Ij— decomposed neutron density pij(r). The filled symbols are the results obtained 
in the continuum HFB calculation while the open symbols are the results obtained in the box- 
discretized HFB calculation. 
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FIG. 4: Neutron pair root-mean-square radius r r .m.s.,«j of the s, p, d and / partial waves of Zr 
isotopes calculated for the ^'-decomposed neutron pair density pij(r). The filled symbols are the 
results obtained in the continuum HFB calculation while the open symbols are the results obtained 
in the box-discretized HFB calculation. 
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FIG. 5: Compositions of different partial waves to the total neutron density, Pij(r)/ p(r), as a 
function of the radial coordinate r for 126 Zr (panel (a)) and 138 Zr (panel (b)). The solid lines 
denote the negative parity states, and the dashed lines the positive parity states. The same plot 
for the neutron pair density, pij(r)/p(r), for 126 Zr (panel (c)) and 138 Zr (panel (d)). 
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FIG. 6: Occupation number densities nij(E) (panels (a) and (c)) and pair number densities nij(E) 
(panels (e) and (g)) of neutron quasiparticle states for different partial waves of 126 Zr, while panels 
(b), (d), (f), and (h) are the same quantities but for 138 Zr. The vertical line indicates the threshold 
energy |A| for the continuum quasiparticle states. 
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FIG. 7: Peak energy E res and the width T of the resonant quasiparticle states around the Fermi 
energy in the P1/2, P3/21 /s/2> an d iV/2 partial waves, plotted in panels (a), (b), (c) and (d) 
respectively, obtained for the Zr isotopes with A = 124 — 138. The filled circles are the resonance 
energy E res , and the vertical bars represent the width T multiplied by a factor of 5. The open 
circles are discretized quasiparticle energies obtained in the box-discretized HFB calculation. 
The dashed line in panels (a) and (b) is the threshold |A| for the continuum quasiparticle states. 
The dash-dotted line in each panel denotes the position Vb of the barrier top of the HF potential 
including the centrifugal potential measured from the Fermi energy. 
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FIG. 8: (a)-(b) Neutron density Airr 2 pij(r) and (c)-(d) neutron pair density 4nr 2 pij(r) of the 
s i/2) Pi/2) P3/2) ^3/2i ^5/2) /5/2; an d /7/2 partial waves in 126 Zr. The total neutron density 
4-Kr 2 p(r) and neutron pair density Airr 2 p(r) are also plotted with the solid line. 
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FIG. 9: (a) (b) Asymptotic exponential constant 2nij of the neutron density pij(r), and (c) (d) kij 
of the neutron pair density pij{r) for the sy 2 , P1/2, P3/2, d 3/2 , d 5 / 2 , / 5 / 2 , and f 7 / 2 partial waves 
in the Zr isotopes. The fitting interval is r = 15 ~ 20 fm. The estimated asymptotic exponential 
constant 2K^ ires and K m i n are also shown with the dashed line and dotted line, respectively, where 
^min = y/4.m\\\/h and Kij tTes = y / 2m(E rcs + \X\)/h 2 are evaluated using the Fermi energy A and 
the resonance energy E rcs shown in Table HI 
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FIG. 10: (a) (b) Neutron density A-7rr 2 pij(r) (solid line), and the truncated neutron density 
47rr 2 p'ij{f) contributed from the low lying quasiparticle states with E < 3 MeV (dashed line) 
(a) for partial wave and (b) for s ly / 2 partial wave in 138 Zr. The contribution 4irr 2 pij(r, E Tes ) of 
the quasiparticle state at the resonance energy E res = 1.057 MeV is also shown (dotted line) for p s / 2 
in panel (a). The asymptotic exponential constant 2ku fitted to the neutron density 4-7rr 2 pij(r), 
2^ fitted to 4-7rr 2 /9^(r), and 2K/j jres calculated with Eq. (|14p for p 3 / 2 partial wave are labeled in 
the corresponding panels, (c) (d) The same as (a) (b) but for the neutron pair densities 4i:r 2 pij(r), 
47rr 2 pij(r) and 4-irr 2 p~ij(r, E vcs ). Note that the unit of the vertical axis is [fm" 1 ] for 4-7rr 2 pij(r), 
47rr 2 p'ij(r), <4itr 2 p~ij(r) and 47rr 2 p\Ar) , but [fm" 1 MeV -1 ] for Anr 2 pij(r, E Tes ) and 4,7rr 2 pij(r, E rcs ). 
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TABLE I: Fermi energy A and the average pairing gaps A uv and A vv for the Zr isotopes. Listed 
are also barrier height of the Hartree-Fock (HF) plus centrifugal potential V ma _ x , HF single-particle 
energies e, resonance energies E Tes and e rcs = E TCS — |A| = E rcs + A, and width T of quasiparticle 
resonances around the Fermi energy for each isotope. All in MeV. 
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-0.446 


-0.427 


-0.415 


-0.401 


-0.380 


-0.336 


-0.232 


-0.126 






0.468 


0.596 


0.656 


0.678 


0.667 


0.628 


0.585 


0.619 




aI 


0.416 


0.532 


0.589 


0.611 


0.605 


0.574 


0.539 


0.577 


Jr i. J Z 




0.289 


0.260 
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0.206 
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0.187 


0.184 
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0.012 


-0.114 


-0.244 


-0.377 
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0.051 
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